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The J/ip yields in heavy ion collisions are expected to be a promising probe of deconfined 
matter, since theoretical models predict that the J/-(/) production could be strongly suppressed 
due to color screening effect in a Quark Gluon Plasma. In addition, one of the interesting 
predictions that has emerged recently is the J/t/j enhancement at RHIC energy via competing 
mechanisms such as recombination. The PHENIX collaboration has measured J/^ production 
via its decay into lepton pairs in Au+Au and Cu+Cu collisions at y^ijvjv ~ 200 GeV. The 
nuclear modification factor Raa for the J/^p is obtained by comparing Au+Au or Cu+Cu 
collisions to p+p collisions. The Raa dependences on centrality, transverse momentum and 
rapidity are presented both at forward rapidity (1.2 < |y| < 2.2) using J/ip tJ-^ fJ-~ and at 
mid rapidity {\y\ < 0.35) using J/tp — > e'^e~ . These preliminary results are compared to cold 
nuclear matter expectations derived from PHENIX d+Au measurements, to results obtained 
by NA50 experiment and to various theoretical models. 



1 3/ip production in heavy ion collisions by the PHENIX experiment 

In heavy ion collisions, the .l/ip production is used as a probe of the produced medium. In 
the case of Quark Gluon Plasma formation, the J/V' yield is expected to be suppressed. The 
NA50 Experiment at the SPS has already measured a suppression in Pb+Pb collisionsl^ this 
suppression is higher than the suppression due to the normal nuclear absorption and is thus called 
anomalous. However, models without QGP formation also describe the observed suppression 
pattern and amplitudel^. At RHIC, the PHENIX experiment aims to observe this effect at a 
10 times higher collision energy. The PHENIX experiment is able to measure J/tp in 2 different 
rapidity regions with magnetic spectrometers: at mid-rapidity in the dielectron decay channel 
{\y\ < 0.35) and at forward rapidity in the dimuon decay channel (1.2 < \y\ < 2.2). Electrons 
are identified in the central arms by their Cerenkov rings and by matching the momentum of 



charged particles reconstructed in drift chambers with the energy deposited in an electromagnetic 
calorimeter. Muons are selected by an absorber and identified by the penetration depth they 
reach in a succession of proportional counters staggered with steel walls. 

At RHIC energy, J/^^'s are mainly produced by gluon fusion. In addition to directly produced 
J/'^'s, measurements also include decays from higher mass Xc and ip' resonances. The i/^ 
production in any A+B ions system is expected to be the production measured in p+p system 
scaled by the number of binary nucleon-nucleon collisions. Thus, results will be shown as 
the nuclear modification factor Rab^ defined as the production measured in A+B collisions 
over the production in p+p times the mean number of binary collisions corresponding to the 
centrality bin used in A+B. However, heavy ion collisions introduce cold nuclear effects such as 
gluon shadowing (modification of the gluon structure function of nucleons in nuclei), transverse 
momentum broadening (the so-called Cronin effect) in the initial state (before the cc formation) 
and nuclear absorption in the final state. Measuring the i/ip yield in d+Au collisions enables 
to quantify these effects. At RHIC energies, results from the d+Au data set^^ indicate a weak 
shadowing and anti-shadowing, which can be described by EKS98 parametrizationl^. The results 
also point out that the nuclear absorption is weak {aabs = 1 to 3 mb) as compared to the one 
measured at the SPS {(Jabs = 4.18 mb). In this presentation, J/-^ results'^ obtained in Cu+Cu 
(3.06 nb~^) and Au+Au (241 /u6~^) collisions are reported and discussed. Data samples from 
d+Au and p+p collisions shown in this paper, have been published by PHENIxl^. 



2 Results and first interpretations 

The nuclear modification factor Rab is presented in figure ^ as a function of centrality (the 
number of nucleon participants Npart) for d+Au, Cu+Cu and Au+Au collisions. In the most 
central Au+Au collisions, the J/V' yield is suppressed by a factor of 3 as compared to the binary 
scaling. This suppression is larger the cold nuclear effect in Au+Au collisions!^, suggesting that 
other mechanisms are involved. 
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Left: Nuclear modification factor Raa as a function of Npart, measured in d+Au (black triangle)^, 
green circles) and Au+Au (gold squares) collisions at toward rapidity 1.2 < \y\ < 2.2 (filled markers) 
and at mid-rapidity \y\ < 0.35 (open markers)^. Black curves correspond to the expected cold nuclear effects in 
Au+Au collisions'^ using EKS98 parametrization and a nuclear absorption cross section aabs = 1 mb (solid line) 
and 3 mb (dashed line). Right: J/ip nuclear modification factor as a function of centrality measured by PHENIX 
(empty circles) and by NA50 (dot)^. The blue dot-dashed line corresponds to the J/t/i expected production 
in NA50((Ta6s — 4.18 mb). The red dashed (solid) line corresponds to the prediction for PHENIX (shadowing 

(EKS98) and nuclear absorption aabs = 3 (1) mb). 

The comparison between NA50 and PHENIX results is an intricate problem. The collision 
energy is 10 times higher at RHIC, leading to a higher temperature which could affect the J/ip 



production differently. Furthermore, the cold nuclear effects may be different: due to a different 
incident energy the nucleus overlap time scale is different between SPS and RHIC, this can mod- 
ify the i /ip nuclear absorption. The different Bjorken x domains also lead to different shadowing 
or anti-shadowing effects. The right figure Q presents the nuclear modification factor Raa as a 
function of the number of participants for NA50 and PHENIX. The suppression amplitudes 
are very similar which is surprising. However, the cold nuclear effects are expected to occur 
differently at the SPS and at RHIC. When using a cross section of 3 mb at RHIC to describe the 
nuclear absorption, the measured over expected ratio amplitude becomes quite similar between 
the 2 experiments, suggesting a possible scenario in which the temperature at RHIC is not high 
enough to reach direct J/V' suppression^. In this case, the measured suppression would be due to 
the sequential suppression of higher mass resonances (xc and ip' decaying into For a more 

detailed comparison, a better understanding of the cold nuclear effects at RHIC is mandatory, 
especially one needs to better constrain the normal nuclear absorption cross section. 

Models without charm quark recombination generally assume that production takes 
place at the early stage of the collisions. In this frame, subsequent interactions in the medium 
(whatever its nature) can only lead to a decrease of the final J/ip yield. At RHIC energy, mul- 
tiple cc pairs are produced in central collisions!^ and quark mobility in a deconfined medium 
allows uncorrelated charm quarks to recombine and to form a J/V'. Such a recombination could 
at least partially compensate 3/ip suppression. However, recombination scenarios imply modi- 
fications of the rapidity and transverse momentum distributions. Recombined J/V' production 
is proportional to the square of the cc pair yields and one expects that recombined production 
will be peaked in the region where the cc pairs are mainly produced: at mid rapidity and at 
low pt according to pQCD predictions. As a consequence, the rapidity distribution and the 
transverse momentum distribution are expected to be narrower in central collisions. The left 
figure 121 presents the rapidity distribution for p+p, Cu+Cu and Au-|-Au collisions for different 
centrality bins: within the error bars, there is no evidence of a modification of the rapidity 
distribution. The right figure [21 presents the (p^) as a function of centrality. The results are in 
between predictions with full and no recombination, suggesting that the measurements consist 
of a mixture of recombined and non recombined J/V'- 

We present here an alternative description of the observed pt broadening using the Cronin 
effect. The broadening is described as a consequence of gluon-ncucleon scattering in the initial 
state. In principle, this effect leads to: 



where {pDpp is the {pD in p-|-p collision, p is the nuclear density, a is the elastic gluon-nucleon 



thickness of nuclear matter of the centrality bin considered. Eq Ogives a reasonable description 
ipl) measurements at lower energy experiments as shown by left figureOlas a function of L. The 
right figure |S1 presents the extrapolation (from p+p and d+Au) of the Cronin effect, showing 
that the pt broadening in Au+Au measured at forward rapidity can be explained with this effect 
alone. This suggests this anomalous J/V' suppression observed in the most central collisions does 
not affect the pt distribution. However, due to error bars, it is difficult to distinguish between 
different models and to make a strong statement about the pt broadening. 




(1) 




"Normalized with p+p collisions studied by NA5lM 

'it is not excluded that direct J/^/) suppression temperature could be reached at RHIC: more statistics would 
allow to use smaller centrality bins and thus to probe a possible direct J/?/) suppression at higher temperature. 
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Figure 2: Left: J/i/) rapidity distribution for p+p, Cu+Cu (red circles) and Au+Au (yellow squares) collisions 
for different centrality bins. Right: (pj) for p+p (black stars), d+Au (blue triangles), Cu+Cu (red circles) 
and Au+Au (yellow squares) as a function of centrality. Filled markers correspond to the forward rapidity 
measurements while open markers to the mid-rapidity measurement. Predictions^^ with full recombination 
scenario for Cu+Cu collisions (red solid line) and for Au+Au (yellow dashed line) and no recombination for 
Cu+Cu (red dotted line) and for Au+Au (yellow dot-dashed line). 
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Figure 3: Left: (p?) as a function of L (fm). Solid bl ue m arkers correspond to Pb+Pb collisions (at 158 G eV 
fixed target) at NA50^| and In+In collisions at NA 60^^ red markers correspond to p+p collisions at NAS^^ 
and p+Cu, p+U, 0+Cu, O+U and S+U at NASS^SI 

(a t 20 GeV fixed target), open blue markers correspond to 

p+Be, p+Al, p+Cu, p+W and p+Pb collisions in 

NA5CPI 

(at 400 GeV fixed target), pink markers correspond to 
p+p et p+d collisions in NA51 and p +Al et p+W collisions in NA50 (at 400 GeV fixed target), orange markers 
are measurements from E866/789/771^kt 800 GeV (fixed target). Right: (pt) as a function of centrality. Red 
solid line (error band in dotted lines) corresponds to the Cronin effect expectation using Eq.^and results from^^. 
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